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Abstract 



Analytic spherically symmetric solutions of the Einstein field equations 
coupled with a perfect fluid and with self-similarities of the zeroth, first and 
second kinds, found recently by Benoit and Coley [Class. Quantum Grav. 15, 
2397 (1998)], are studied, and found that some of them represent gravitational 
collapse. When the solutions have self-similarity of the first (homothetic) kind, 
some of the solutions may represent critical collapse but in the sense that now 
the "critical" solution separates the collapse that forms black holes from the 
collapse that forms naked singularities. The formation of such black holes 
always starts with a mass gap, although the "critical" solution has homothetic 
self-similarity. The solutions with self-similarity of the zeroth and second 
kinds seem irrelevant to critical collapse. Yet, it is also found that the de 
Sitter solution is a particular case of the solutions with self-similarity of the 
zeroth kind, and that the Schwarzschild solution is a particular case of the 
solutions with self-similarity of the second kind with the index a = 3/2. 
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Typeset using REVT^ 



I. INTRODUCTION 



Gravitational collapse of a realistic body has been one of the most important and thorny 
subjects in General Relativity (GR) since the early times of GR Quite recently, thanks to 
Choptuik's numerical discovery of critical phenomena in the threshold of black hole formation 
0], the subject has attracted further attention. As a matter of fact, it is so attractive that 
Critical Phenomena in Gravitational Collapse has been already a very established sub-area 
in GR [0,0. From all the work done so far, the following seems clear: (a) The critical solution 
and the two dimensionless constants A and 7 are universal only with respect to the same 
matter field, and usually are matter-dependent, (b) The universality of the critical solution 
and the exponent 7 now are well understood in terms of perturbations , while the physical 
origin of A still remains somewhat of a mystery. The former is closely related to the fact 
that the critical solution has only one unstable mode. This property now is considered as 
the main criterion for a solution to be critical, (c) The critical solutions can have discrete 
self-similarity (DSS) or homothetic self-similarity (HSS) [|, or none of them, depending on 
the matter fields and regions of the initial data spaces. So far, in all the cases where the 
critical solution either has DSS or HSS, the formation of black holes always turns on with 
zero mass, the so-called Type II collapse, while in the cases in which the critical solution has 
neither DSS nor HSS, the formation always turns on with a mass gap, the so-called Type I 
collapse 

In the Type II collapse, it is the usual belief that the fact that black hole starts to form 
with an infinitesimal mass is closely related to the fact that the problem concerned is of 
scale invariance, for example, the Einstein equations coupled with a massless scalar field. 
When the scalar field is massive, the corresponding field equations are scale-invariant only 



In the literature, homothetic self-similarity has been also called continuous self-similarity. How- 
ever, in order to distinguish it from the self-similarity of the other kinds, in this paper we shall 
refer it as homothetic self-similarity, or self-similarity of the first kind. 
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asymptotically . For a perfect fluid with the equation of state p = kp, the corresponding 
Einstein field equations are also of scale invariance. As a result, in all these cases critical 
phenomena of Type II collapse were found, and in the case of the scalar field the critical 
solution has DSS, while in the case of perfect fiuid, the critical solutions have HSS 

It is known that homothetic self-similarity is a particular case of kinematic self-similarity 
P|. In fact, the latter consists of three kinds, the zeroth kind, the first (homothetic) kind, 
and the second kind. Thus, a natural question is: Can critical solutions have self-similarity 
of the other kinds? 

In this paper, we shall study this problem for the gravitational collapse of perfect fiuid 
with kinematic self-similarity. As a matter of fact, several classes of such analytic solutions 
to the Einstein field equations are already known p. So, here we shall study these solutions 
in some details and pay particular attention on critical solutions. Finding critical solutions 
usually consists of two steps, one is first to find a generic family (or families) of solutions, 
characterized, say, by a parameter p, such that when p > p* the collapse forms black holes, 
and when p < p* it does not. Once such solutions are found, one needs to make perturbations 
of the solution p = p* and to study the spectrum of their modes. If the solution has only 
one unstable mode, then by definition this solution is a critical solution, and the exponent 
7 is given by 

7 = r^, (1.1) 

l^^il 

where cxi is the unstable mode . In this paper, we shall consider only the first part of the 
problem, and leave the study of perturbations to another occasion. Specifically, the paper is 
organized as follows: In Sec. II we shall give a brief introduction to kinematic self-similarity, 
and in Sec. Ill we shall study the Benoit-Coley (BC) solutions with self-similarity of the 
zeroth kind, while in Sec. IV the BC solutions with self-similarity of the first and second 
kinds will be studied. The paper is closed with Sec. V, in which our main conclusions are 
presented. An appendix is also included, where the Einstein field equations are written in 
terms of self-similar variables. 
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II. SPHERICALLY SYMMETRIC SPACETIMES WITH KINEMATIC 

SELF-SIMILARITY 



Self-similarity refers to the fact that the spatial distribution of the characteristics of 
motion remains similar to itself at all times in which all dimensional constant parameters 
entering the initial and boundary conditions vanish or become infinite [0. Such solutions 
describe the "intermediate asymptotic" behavior of solutions in the region where a solution 
no longer depends on the details of the initial and/or boundary conditions. 

Cases in which the form of the self-similar asymptotes can be obtained from dimensional 
considerations are referred to as self-similarity of the first (homothetic) kind [|10|- Solutions 



of the first kind were first studied by Cahill and Taub in GR for a perfect fluid [|rT|. They 
showed that the existence of self-similarity (of the first kind) could be formulated invariantly 
in terms of a homothetic Killing vector, S,^, which satisfies the conformal Killing equation, 

^^9i^u = 2gf,u, (2.1) 

where C denotes Lie differentiation along From the above it can be shown that 

A^Mv = 0. (2.2) 

For a perfect fluid with the energy-momentum tensor (EMT) given by 

Tf,u = {p + p)Uf,Uy - pg^y, (2.3) 

it can be shown that it is consistent with Eq. (|2.2| ) if we require 

C^u^" = -u^, C^p = -2p, C^p = -2p. (2.4) 

Hence, in this case "geometric" self-similarity and "physical" similarity coincide, although 
this does not need to be so in more general cases |8|JT^ . Applying the above to the spacetimes 
with spherical symmetry. 



where dfl"^ = dO"^ + sin^ Odip^ , Cahill and Taub found that the condition (|2.1|) requires 

$(t,r) = $(0, vI/(t,r) = vl/(^), s{t,r) = S{0. (2-6) 

where 

r 

The corresponding homothetic Kilhng vector is given by 

d d d 

^ dx^^ dt ~^ dr ^ 

Note that in writing the metric ( p.5|) we had multiphed a factor rf to the usual spherical 
metric, so that the metric coefficients $, S and the coordinates t, r, 9 and now are all 
dimensionless, where we assume that ri has the dimension of length. It is found that this 
choice will simplify the dimensional analysis to be given below. The corresponding Einstein 
tensor and Einstein field equations are given in terms of both t, r and x, t in the Appendix, 
where x and r are the self-similar variables that are functions of t and r. Their explicit 
definitions in each case are given in the Appendix. 

The existence of self- similarity of the first kind is closely related to the conservation laws 
and to the invariance of the problem with respect to the group of similarity transformations 
of quantities with independent dimensions, in which case a certain regularity of the limiting 
process in passing from original non-self-similar regime to the self-similar regime is assumed 
implicitly. However, in general such a passage does not need to be regular. Consequently, 
the expressions for the self-similar variables are not determined from dimensional analysis. 
Such solutions are then called self-similar solutions of the second kind. A characteristic of 
these solutions is that they contain dimensional constants that are not determined from the 
conservation laws [0. Using these arguments to a perfect fiuid ( |2.3| ), Carter and Henriksen 



g] gave the notion of kinematic self- similarity with its properties, 

r^V = A^'' = (2-9) 

where h^i, is the project operator, defined by 



h^,u = gflu - u^Uu, (2.10) 

and a is an arbitrary dimensionless constant. When a = 1, it can be shown that the kine- 
matic self-similarity reduces to the self-similarity of the first kind (homothetic self-similarity). 
When a ^ 1, Carter and Henriksen argued that this would be a natural relativistic coun- 
terpart of self-similarity of the second kind {a ^ 1), and of the zeroth kind {a = 0), in 
Newtonian Mechanics. 

Applying the above to the spherical case, Carter and Henriksen found that the metric 
coefficients $, \E' and S should also take the form of Eq.( p.6|) but with the self-similar 
variable ^ and conformal vector now being given, respectively, by 



for the second kind, and 



E^4r' « = = (2.12) 



for the zeroth kind. Comparing Eg. ( f^.8[ ) with Eq.( p^.ll[ ) we find that the self-similarity of 



the first kind can be considered as a particular case of the one of the second kind. In this 
paper we shall do so, although, as we mentioned above, the physics is quite different in the 
two cases. In particular, when the coordinates t and r are rescaled, t' = ct and r' = cr, 
where c is a constant, C, is unchanged only for the homothetic case a = 1. 

III. ANALYTIC SOLUTIONS OF PERFECT FLUID WITH SELF-SIMILARITY 
OF THE ZEROTH KIND AND THEIR PHYSICAL INTERPRETATIONS 

The solutions to be studied in this section are those given by Eqs.(2.54) - (2.57) in 0. 
Note that the expression for the function S{x) given there is not correct. As a matter of 
fact, setting $ = 0, we find that Eq. ( |A.27D yields 

e'" = {l + yfS', (3.1) 
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while Eg. ( |A.2^) is satisfied automatically. Submitted Eq.( p.l|) into Eq. ([A.26|) , we obtain 



y,xx + ^yy,x = o, (3.2) 

which has the first integral 

2y,, + 3y^+Po = 0, (3.3) 

where x = In(^), and po is a dimensionless constant, in contrast to the claim given in [§. 
Then, it can be shown that the corresponding perfect fluid is given by 

P = ^2fr, — T' -P = — ' = '"I'^M' 3.4 
rf[l + y) rf ^ 

where y{x) = S^^l S. Depending on the sign of po; Eq.( p.3|) has physically different solutions. 
In the following let us consider them separately. 

A. Case po = 
In this case, it can be shown that Eq. (|3.3|) has the solution, 

where 5*0 and Xq are integration constants. Without loss of generahty, we can set 5*0 equal 
to one by a conformal transformation. In the following we shall assume that this is always 
done whenever it is applicable. Substituting the above expressions into Eq. ( |3.4|) , we find 
that 

_ 4 

^" r2(x + Xo) [3(x + Xo) + 2]' ^ " °' ^^'^^ 

which shows that in this case the solutions represent a dust fluid. Thus, these solutions 
must belong to the general Tolmann-Bondi class . 



^In this paper the units will be chosen such that the Einstein coupling constant k [= SvrG/c^] = 1. 
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FIG. 1. The spacetime described by the solutions of Eq.(|3.5|) in the text in the (t,r)-plane. 
It is singular on the hypersurfaces x + xq = and x + xq = —2/3, which divide the whole 
spacetime into three disconnected regions: / = {x^ : x + xq > 0} , // = {x^ : —2/3 < x + xq < 0}, 
and 111 = {x'^ : X + xo < -2/3}. 

From Eq. ( |3.(j| ) we can see that the spacetime is singular on the hypersurfaces x + xq = 
and x + Xq = —2/3. These two hypersurfaces divide the whole spacetime into three regions, 
I, II, and III, where / = {x^ : x + Xq > 0} , // = {x^ : -2/3 < x + Xq < 0}, and /// = 
{x^ : X + Xq < — 2/3} [See Fig. 1]. In region II, the energy density of the fluid p is negative, 
and the physics of the spacetime in this region is not clear. In region III, it is non-negative 
and the singularity located on the hypersurface x + xq = —2/3 is naked, and the spacetime 
in this region can be considered as representing an inhomogeneous cosmological model. In 
region I, to study the nature of the singularity located on the hypersurface x + Xq = 0, let 
us first calculate the gradient of the geometric radius, R = rS{x), of the two sphere, 

= 1 ^2/3 W - + ^o)'/'} • (3.7) 
9(x + xoj ' 

The formation of apparent horizons are indicated by the vanishing of the gradient. Thus, 
setting the right-hand side of Eq.( |3.7|) to zero, we obtain 
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X + Xo = ) ' , {RaRf39''^ = O) . (3.8) 

Clearly, for any given r, we always have x + xq > on the apparent horizon. Hence, in 
the present case the formation of the spacetime singularity on the hypersurface x + xq = 
always follows the formation of the apparent horizon, or in other words, the singularity is 
always covered by the apparent horizon. These solutions can be considered as representing 
the formation of black holes due to the gravitational collapse of the fluid, starting at the 
moment t = — oo. Defining the mass function m{t,r) as fl^. 



m(t,r) = |(l + i?,„i?,^^7°^), (3.9) 

we find that m{t, r) = 2r^/9. On the apparent horizon, we have 

MBH{t)=m{t,rAH{t)) = (3.10) 

where rAuif) is a solution of Eg. (|3.8|). The quantity MBnif) can be considered as the contri- 
bution of the collapsing perfect fluid to the total mass of such formed black holes [0], which 
in the present case goes to infinity as rAH{t) +oo, as can be seen from Eq.( p.8|) . There- 
fore, now the collapse of the dust fluid always forms black holes with infinitely large mass. 
To remedy this shortage, one may cut the spacetime along a non-spacelike hypersurface, 
say, r = ro(t), and then join the region, r < ro(t), to an asymptotically flat region |15|. By 



this way, we can see that the resultant model will represent the collapse of a dust ball with 
its radius r = ro(t) [cf. Fig. 2]. From the moment t = tf on, the ball collapses completely 
inside the apparent horizon, and the contribution of the collapsing ball to the total mass of 
such a formed black hole is given by 

MBH=MBH{tf), (3.11) 

where tj is a solution of the equation, 

ro{tf) =rAH{tf). (3.12) 
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In the present case, since the fluid is co-moving with the coordinates, without loss of gener- 
ahty, we can choose the joining hypersurface as ro(t) = = Const. Then, from Eq.( p.ll[ ) 
we flnd that Mbh{= 2rJ/9) is always finite and different from zero for any given non-zero 
tq. This is different from the gravitational collapse of dust fluid with self-similarity of the 
flrst kind studied in |TB|, where it was shown that for any given non-zero Tq, black holes 
with inflnitesimal mass can be formed by properly choosing a parameter that characterizes 
the strength of the collapse. 



R = 




FIG. 2. The corresponding Penrose diagram of the spacetime described by Eq.( |3.5| ), after it 
is first cut along the hypersurface r = tq and then joined with an asymptotically flat region. The 
point P represents the moment when rg = rAH{ij)) where the ball of the fluid collapses completely 
inside the apparent horizon. 



B. Case po > 

In this case it can be shown that Eq.( p.3|) has the solution. 
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y{x) = — atanA(x), S{x) = cos^^^ A{x), 



(3.13) 



where Xq is another integration constant, and 



3a 



A{x) = — {x + Xq) , a 



3 



1/2 



(3.14) 



Then, from Eg. (|3.4| ) we find that 



X + Xq 



_2_ 

3a 

27r 
3a 



-1/1 



3a^tany4(a:;) [tany4(a;) — a] Po ^ 

^ rj[l — ata.nA{x)] ' ^ rj 

from which we can see that the solutions are singular on the hypersurfaces, 

tan" 
(n + i) 

except for the case a = 1, where we have 

p = -■^tan|^(x + xo) }> , (a = l). 



(3.15) 



(3.16) 



(3.17) 



where n is an integer. Clearly, in the latter case the spacetime is singular only on the 
hypersurfaces 



27r / 1\ , 
a; + a:o = Y U+ - I , (a = 1). 



(3.18) 



However, in either of the two cases, the spacetime is singular on various hypersurfaces, and 



the energy conditions, weak, strong and dominant IT^, hold only in certain regions. The 
physics of these singularities are not clear, and the solutions may have physical applications 
only in certain regions. In particular, they cannot be interpreted as representing gravita- 
tional collapse of the fiuid. 



C. Case Po <0 

In this case, the solutions can be further classified into three different cases, according 
to a) y"^ > a^, j3) y"^ = a^, and 7) y"^ < . 

a) Case > c?: In this case it can be show that Eq. (|3.3| ) has the solution, 
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smn A[xj 



(3.19) 



while Eqs. (|3.4|) yields, 



P 

P 



3a^ coshA(x) [coshy4(a:) + asinliy4(x)] 
rf sinh A{x) [sinh A{x) + a cosh A{x)] ' 
\po\ 3a2 



^2 



(3.20) 



where A{x) and a are still given by Eq.( p.l4D . The three energy conditions now require 
p > 9a^. Then, from Eq. (|3.2(]| ) we can see that this condition holds only in the region 



Inr + xo — xi < t < Inr + xq, 



(3.21) 



for any given a, or in the region 



t<lnr + Xo + X2, (a < 1), 



(3.22) 



for a < 1, where 



Xi 



X2 



3a \ 1 + a I 



3a 



1 + a 



> 0. 



(3.23) 



In the region defined by Eq.( p.2lD , the spacetime is limited by the curvature singularity 
located at a; + xq = in one side, and by the hypersurface x + xq = xi in the other side, 
across the latter the energy conditions do not hold. The solutions in this region seem not to 
have much physics. In the region defined by Eq.( p.22| ), the spacetime may be considered as 
representing a cosmological model. It is interesting to note that the spacetime in this region 
is free of singularities and asymptotically fiat as t ^ +oo. However, it is not geodesically 
complete and needs to be extended beyond the hypersurface x + Xq = —X2. A "natural" 
extension would be the one simply given by the above solutions (|3.19| ). This extension will 
be valid until the hypersurface x + Xo = —x^, where 



X3 = -—In 
3a 



1 — a 



1 + a 



(3.24) 
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on which the spacetime is singular. Obviously, the fluid in this extended region do not 
satisfy all the three energy conditions. 

j3) Case y"^ = In this case, it can be shown that Eq.( |3.3| ) has the solution. 



3^2 

y{x) = ±a, S{x) = e"^"^, p = -p = 



(3.25) 



Introducing a new radial coordinate r via the relation, r = r^^", we flnd that the corre- 
sponding metric can be written in the form. 



which is the de Sitter solution [|T7 . 

7) Case y"^ < a?: In this case, we flnd that 



cosh A{x) 



X 



3a^ sinh A{x) [sinh A{x) + a cosh A{x)] 
rfcoshyl(x) [coshA(x) +asinhA(x)] ' 
3a2 



P 



It can be shown that now the energy conditions hold only in the region. 



In r + Xo + ^3 < t < In r + Xo + 3:4, (a > 1) 



(3.26) 



(3.27) 



(3.28) 



for a > 1, where X3 is given by Eg. ( p. 24] ) and X4 is given by 

a — 1 



Xi 



3a 



(3.29) 



2 + VsTo^ 

When a < 1, there does not exist any region in which the three energy conditions hold. The 
spacetime is singular on the hypersurface x + Xq = — X3. The physics of the spacetime in 
this case is not clear (if there is any). 
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IV. ANALYTIC SOLUTIONS OF PERFECT FLUID WITH SELF- SIMILARITY 
OF THE FIRST AND SECOND KINDS AND THEIR PHYSICAL 

INTERPRETATIONS 



The solutions to be studied in this section are those given by Eqs.(2.27) - (2.31) in [Q, 
for which we have $ = 0. From Eq.( A.45|) we find that the function \1/ takes the same form 



as that given by Eq.(|3.1|) in terms of y, while Eq. (|A.4^ yields 



y,xx + (3y + = 0, (4.1) 

which allows the first integral, 

2y,.. + + 2ay + a V = o, (4.2) 

where po is also a dimensionless constant, in contrast to what claimed in 0. It can be shown 
that in the present case Eq.( [A.39D is satisfied automatically, too. Then, the corresponding 



perfect fiuid is given by 

y [(3 - 2a)y - a'po] 

When a = 1, the corresponding solutions have self-similarity of the first kind, otherwise, 
they have the second kind. 

Note that the solutions of Eq.( [4.2D for the function y{x) given in are not correct. 
Thus, in the following we shall first derive the correct expressions for y{x) and S{x), and 
then study the physics of the solutions. Depending on the value of po? the solutions of 
Eq.( [4.2|) can be divided into several classes. In the following let us consider them one by 
one. 

A. Case po = 
When pq = 0, Eq.(ff]2|) has the solution. 
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3 [e"{^+^o) _ 1] ' 



S{x) 



^a{x+xo) _ ]^ 



2/3 



(4.4) 



where Xq is another integration constant. The corresponding energy density of the fluid is 
given by 



P 



4(3 - 2a) r . 



9rf 



^a{x+xo) _ 



^a{x+xo) _|_ 



2a - 3 



(4.5) 



Since p = in the present case, the above solutions must also belong to the general Tolmann- 
Bondi solutions |T^. 

a) Case < a < 1: In this subcase, Eq.( [4.5| ) shows that the spacetime is singular on 
the hypersurfaces. 



a) t = 0, b) X + Xq = 0, c) X + Xq = — X5, 



(4.6) 



where x^ is defined as 



x^ = In 

a 



3 -2a 



> 0. 



(4.7) 



From the above we can show that p is non-negative only in the region x + xq < —x^ or 
in the region x + Xo > 0. In the region x + xq > 0, the spacetime is singular on the two 
hypersurfaces x + xq = and t = 0. While the physical meaning of the spacetime in the 
region x + xq > 0, and t < is not clear, the spacetime in the region t > can be considered 
as representing a cosmological model with its initial singularity at t = 0. The spacetime 
in the region x + Xq < —x^ can be considered as representing the gravitational collapse of 
the perfect fluid. To study the nature of the spacetime singularity at x + xo = —x^, let us 
consider the quantity, 

4Q;^e^f^"^~"^^~"^°' 



from which we find that 

{-tAH)~ -- 



9(_t)2(a-l)/Q ^Q-a{x+xo) _ 
2^g{l-a)x-axo 



2/3 



(4.8) 



3 ^^-a(x+xo) _ 1] 
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1/3, [R,aR,f3g"''^ = 



(4.9) 



FIG. 3. The spacetime described by the solutions of Eq.(^!J) for < a < 1 in the (t,x)-plane. 
It is singular on the hypersurfaces o) t = 0, b) x = —xq and c) x = — (xq + 2:5). At the moment 
t = tc the apparent horizon crosses the singular hypersurface x = — (xq + X5), and asymptotically 
approaches to the point (t,x) = (0, —xq). 

It can be shown that this hypersurface will cross the singular hypersurface x + xq = —x^ 
at the moment t = tc, where tc is given by 

i-tc)'-^ = ^ [e""^ - 1]'/^ e(^-")"^+"°. (4.10) 

This can be seen clearly in the (t, a;)-plane, as illustrated by Fig. 3. The corresponding 
Penrose diagram is given by Fig. 4, from which we can see that the spacetime singularity 
formed at x + Xq = — X5 is covered by the apparent horizon at the beginning {t < tc)- As 
the fluid continues to collapse, the apparent horizon starts to form after the formation of 
the spacetime singularity, hence it becomes naked when t > tc. 



17 




FIG. 4. The Penrose diagram for the spacetime described by the solutions of Eq.([4.4|) for 
< a < 1. The singularity formed on the hypersurface x + xq = —x^ is covered by the apparent 
horizon when t < tc and becomes naked when t > tc, where tc is defined by Eq.( 4.1(]| ). 



b) Case a = 1: In this subcase, it can be shown that the apparent horizon is given by 

x + xo = -xe, (RaRpg''^ = O) , (4.11) 




where 

X6 = ln|l + ^e-=^^o^. (4.12) 

Since 

/27 + 8e~^^'"\ 

xe-x, = \ni j = {=0, Xo=p% (4.13) 



where p* = — [ln(27/4)]/3, we find that, when xq < p*, the apparent horizon always forms 
before the formation of the spacetime singularity at x + xq = — X5, that is, the collapse now 
always forms black holes. When xq = p*, the apparent horizon and the spacetime singularity 
are formed on the same hypersurface, i.e., now the singularity is marginally naked. When 
Xo > p*, the apparent horizon always forms after the formation of the spacetime singularity, 
or in other words, now the collapse always forms naked singularities. The contribution of 
the collapsing fluid to the total mass of such formed black holes is given by 



2p-2a:o 

MBH{t) = ——TAHit), (4.14) 

y 



where rAuit) is a solution of Eg . (14.111 ). Thus, as r^_H-(t) +00, we find Mbh +00. 



Similar to the case discussed in the last section, to obtain a black hole with finite mass, 
we can cut the spacetime along the hypersurface r = tq = Const, and then join the region 
r < Tq with an asymptotically fiat region. By this way, the resulting model will represent 
gravitational collapse of a ball with its comoving radius r^. At the moment t = tj, where tf 
is a solution of the equation TAuitf) = "'^O; the ball collapses completely inside the horizon, 
and its contribution to the total mass of such formed black holes is given by 

2p-2xo 

Mbh = —^^0, (4.15) 

which is always finite and non-zero for any given non-zero tq. It is interesting to note that in 
the present case the solutions may represent critical phenomena. To have a definite answer 
to this problem, we need to study the spectrum of perturbations of the "critical" solution 
and show that it has only one unstable mode. This is currently under our investigation. 

It is very interesting to note that in this case the black holes start to form with a mass 
gap, although the Einstein field equations are of scale invariance, and the spacetime has 
self-similarity of the first kind [^Q. Thus, the solutions studied in the present case show 
clearly that even the solutions have homothetic self-similarity, the formation of black holes 
does not necessarily start with an infinitesimal mass. 
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t 




FIG. 5. The spacetime described by the solutions of Eq.(4.4) for 1 < a < 3/2 in the 
(t, a;)-plane. It is singular on the hypersurfaces a) t = 0, b) x = —xq and c) x = —{xq + X5). 
At the moment t = tc the apparent horizon crosses the singular hypersurface x = — (xq + 2:5), and 
asymptotically approaches to the one x = —xq. 



c) Case 1 < a < |: In this subcase, it can be shown that the spacetime is also singular 



on the hypersurfaces given by Eq.(4^) and the apparent horizon is given by Eq. (|4.9| ). In 
the (t, x)-plane, it is given by Fig. 5, from which we can see that it also crosses the singular 
hypersurface x + Xq = — X5 once, but in contrast to the subcase < a < 1, now the 
singularity initially is naked and becomes covered by the apparent horizon after the moment 
t = tc, as shown by Fig. 6. 



d) Case a = |: From Eq.(I]5) we find that p = = p. That is, in this subcase the 



spacetime is vacuum. The metric coefficients are given by 



y{x) 



S{x) 



-{x+xo) 



g3{x'+xo)/2 _ ^ 



2/3 



(4.16) 



Defining a new radial coordinate f via the relations. 



f ^ 1 g3xo/2^3/2 



(4.17) 



we find that the metric can be written as 
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where 



ds = dr 




- r 
2^ 



4/3 



d^^ 



(4.18) 



r 



-t, Tg = e 



-3x0 



(4.19) 



This is exactly the Schwarzschild solution written in the Lemaitre coordinates |]T^, with Vg 
being the Schwarzschild radius. From this case we can see that the parameter is related 
to the total mass of the Schwarzschild black hole. 




FIG. 6. The Penrose diagram for the spacetime described by the solutions of Eq.([4.4|) for 
1 < a < 3/2. The singularity formed on the hypersurface x + xq = —x^ is naked for t < and 
covered by the apparent horizon when t > tc, where tc is defined by Eq.( [4.10| ). 



e) Case a > |: In this subcase the energy density of the fluid takes the form 



4(2a - 3 ) 
9rf 



1-e 



a{x+X()) 



2a - 3 



(4.20) 



Thus, p > requires 
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x + xo<0, (p>0). 



(4.21) 



















FIG. 7. The spacetime described by the solutions of Eq.(^) for a > 3/2 in the {t, x)-plane. It 
is singular on the hypersurface x = —xq and the energy density of the fluid is non-negative only in 
the region x < —xq. 



R = 




FIG. 8. The Penrose diagram for the spacetime described by the solutions of Eq.([4.4|) for 
a > 3/2. The apparent horizon now always forms before the formation of the spacetime singularity 
at X + Xq = 0. 
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In this region the spacetime is singular only on the hypersurface x + Xq = 0. On the other 
hand, the apparent horizon in the present case is still given by Eq.( ^.9| ). In the (t, x)-plane, 
this hypersurface is shown by Fig. 7. The corresponding Penrose diagram is given by Fig. 8, 
from which we can see that now the apparent horizon always forms before the formation of 
the spacetime singularity, that is, the solutions now represent the formation of black holes. 
The contribution of the fluid to the total mass of such formed black holes is given by 



MBHit) 



(4.22) 



where rAH and xah are the solution of Eq. (|4.9|) . When the spacetime is first cut along the 
hypersurface r = tq and then joined with an asymptotically flat region, and the contribution 
of the fluid to the total mass of such formed black holes is given by 

-|2/3 



M 



BH 



2 



-a{xf+xo) 



(4.23) 



where Xf = ln[ro/(— t/)^/""], and t/ denotes the moment when the ball collapses completely 
inside the apparent horizon, given by tq = rAuiif)- One can show that this mass is also 
finite and non-zero for any given non-zero tq. 

B. Case < po < 5 
In this case, it can be shown that the corresponding solutions are given by 



-(X + Xoj 



a 



S{x) = e"""/^ cosh'/^ 



X + Xq) 



(4.24) 



where /? = a|l — 3po|^^^- Then, the energy density of the fluid is given by 



P 



(3{2a - 3) • 


(1 -3po)'^/^ - tanh 


f (x + Xo) 


} 


Sa'^rjt'^ |tanh 


|(x + Xo) 


.A} 



X < tanh 



+ B 



(4.25) 



where 
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A 



a 



n Q |i/2' 
a |1 — 3po 



B 



3apo — (2a — 3) 

n72 



(4.26) 



^01 ■ |l-3po| ^ (2a-3) 

From the above equation it can be shown that, when < a < ai, we have A > 1; when 
O-i < a < we have —1<A< +1; and when a > a2, we have A < —1, where 

3 3 



l + (l-3po)i/^ 



a2 



When < a < 3/2, we have B < -{1 - 3pc 



1 



v-l/2 



3po)i/2' 



(4.27) 



< — 1; when 3/2 < a < «i, we have 



B > +1; when ai < a < a2, we have —1<B< +1; and when a > 0^2, we have B < —1 
[cf. Fig. 9]. Thus, from Eq. (|425| ) we find that p is always non-negative at any given point 
of the spacetime when < a < ai or when a > 02- It is singular on the hypersurface t = 0. 
This hypersurface divides the spacetime into two disconnected regions, t > and t < 0. 
In the region t > 0, the spacetime can be considered as representing a cosmological model 
with its initial singularity at t = 0. In the region t < 0, the spacetime can be interpreted 
as representing the gravitational collapse of the perfect fiuid. In this region, an apparent 
horizon is formed on the hypersurface, given by 



cx — I 

[-tAH)~ 



ae 



{3-a)x/3 



3 cosh^/^ 



|(x + Xo) 



X < cosh 



-{X + Xo) 



:i-3p,y/'^,{R^Rf,g'^^ = 0). 

(4.28) 
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FIG. 9. The curves of the functions A{a) and B{a), defined by Eg. (4.26) versus a 



It is not difficult to see that this hypersurface is formed always before the formation of 
the spacetime singularity at t = 0. Thus in the present case the collapse always forms black 
holes, and the contribution of the fluid to the black hole mass is given by 



-^{xah + Xo) 



(4.29) 



where rAuif) and xah are the solutions of Eq.( |4.28| ). Similar to the cases discussed above, 
the mass of such formed black holes now becomes also infinitely large as VAHit) +oo. 
Thus, in this case we also need to cut the spacetime along the hypersurface r = tq and then 
join the region r < vq to an asymptotically fiat region. Once this is done, it is not difficult 
to see that 



2^ + ^oJ 



-axf /3 



(4.30) 



where Xf is given by Xf = ln[ro(— t/) ^/"], and tf denotes the moment when the ball of 
perfect fiuid collapses completely inside the apparent horizon, which is given by Eq.( [4.28|) 
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with TAnitf) = To- Clearly, for any given non-zero tq, Mbh is non-zero. That is, in the 
present case the black holes start to form with a mass gap, too. 




FIG. 10. The spacetime described by the solutions of Eg. ( 4. 24 ) for ai < a < 02 in the 
{t, x)-plane. p is non-negative only in the region x > — (xq + X7) or in the region x < —{xq + xg). 
It is singular on the hypersurface x = — (xq + X7). 



When ai < a < 0:2, p is non-negative only in the region x + xq > —x-j or in the region 
X + xo < —xs, where 



2-1 2 
a;7 = — tanh (A), = -5 tanh (B). 

P P 



(4.31) 



Since B > A for ai < a < 0^2, we find Xs > Xj, where equality holds only when a = ai, 
or a = a2- The spacetime is singular on the hypersurface t = and x + xq = —X7. Once 
again, the region t > can be considered as representing a cosmological model. In the 
region x + Xq < — X7, on the other hand, we find that 

— 00, X + Xq = —Xj, 
< 0, —Xs < X + Xq < —Xy, 

(4.32) 

= 0, X + Xq = —Xs, 
. > 0, X -|- Xo < — Xg. 

In this region, the apparent horizon is still given by Eq.( [4.28D , from which we find that, as 
X ±00, we have (— j(:)("^i)/° +00. In the (t,x)-plane, it is given by a curve that crosses 



P 
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both the hypersurfaces x + Xq = —Xj and x + Xq = —Xs [cf. Fig. 10]. The corresponding 
Penrose diagram is given by Fig. 11, the physics of which is unclear. 




FIG. 11. The Penrose diagram for the spacetime described by the solutions of Eq.( [4.24| ) for 
0(1 < a < 02- In the region — xy < x + xq < —xs, the energy density of the fluid becomes negative. 

C. Case po = h 



In this case, it can be shown that Eqs. (|4.2| ) and (|4.3| ) have the following solutions, 

y = Tr^-^^ s = {x + xor'e-^i\ 

3(x + Xo) 3 

1 

P 



3a2rft2(x + xo) [(3 - a){x + xq) + 2] 
X |a^(3 - a){x + xof 



-6a(2 -a){x + xq) + 4(3 - 2a)} 
1 

3^t2 



P=^2- (4-33) 



A) Case < a < 1: In this subcase from Eq.( [4.33D we can see that the spacetime is 
singular on the hypersurfaces 
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a)t = 0, b)x + XQ = 0, c)x + Xo = -- . (4.34) 

3 — a 

It is not difficult to show that now the singularity at x + xq = —2/ (3 — a) is first formed 
and the one at t = is last formed. The spacetime in the region t > may be considered as 
representing cosmological model with its initial singularity of the spacetime at t = 0. The 
region x + xq < —2/ (3 — a) can be considered as representing the gravitational collapse of 
the perfect fluid. To study the nature of this singularity, let us consider the formation of 
apparent horizons, given by 

In the {t, x)-plane, this curve is similar to that given in Fig. 3, if x^ is replaced by 2/(3 — a). 
Thus, in this case the singularity at x + xq = —2/(3 — a) is covered upto the moment t = tc, 
where now is given by 

(-g^ = 1 (i^)^'%V3. ,4.36) 

After this moment, the singularity becomes naked. The corresponding Penrose diagram is 
given by Fig. 4. 

B) Case a = 1: As we mentioned previously, when a = 1, the corresponding solutions 
have the self-similarity of the first kind. Setting a = 1 in the above expressions, we find 
that 

^ {X + Xq)^ - 3{X + Xq) + 2 ^ J_ 

^ 3rlt^ix + Xo)[{x + Xo) + iy ^ 3r?t2- ^ 

Clearly, the spacetime is also singular on the hypersurfaces given by Eg. ( [4.34] ). Since now 
we have = r/{—t), we can see that these singular hypersurfaces are straight lines in 
the {t, r)-plane. The spacetime in the region t > may be interpreted as representing 
cosmological model with an initial spacetime singularity at t = 0. The region x + xo < —1 
may be considered as representing the gravitational collapse of the perfect fluid starting at 
t = — oo. To study the nature of the singularity located on the hypersurface x + Xq = — 1, 
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let us, following the analysis given above, first consider the formation of apparent horizon 
in this region, 

4x /3 

where 

Yi{x) = [2-{x + Xo)f , ^2(2;) = 9{x + a;o)2/'e-^^/3_ (4 39) 

It is easy to show that Y2{x) has one minimal at x + xq = and one maximal at x + xq — 1/2. 
When X ^ —00 it diverges exponentially, and when x — > +00 it goes to zero exponentially. 
On the other hand, Yi{x) is a parabola with its minimum located at x + = 2 [cf. Fig. 
12]. Thus, in general the equation Yi{x) — Y2{x) has three real roots, say, xg, xiq and 
Xii. Without loss of generality, we assume that Xu > Xiq > Xg. Then, we can see that 
X — xg represents the outmost trapped surface, i.e., the apparent horizon. Introducing a 
new parameter D via the relation, 

D=-{xg + xo + l), (4.40) 

we find that the hypersurface x = Xg can be written as 

x + xo^-{l + D), {RaRpg"^ = 0) . (4.41) 

Thus, when D > 0, the apparent horizon always forms before the formation of the spacetime 
singularity at x + xq = —1. That is, in this case the gravitational collapse of the perfect 
fluid always forms black holes. When D < 0, the apparent horizon always forms after the 
formation of the spacetime singularity at x + Xq — —1, namely, now the collapse always 
forms naked singularities. When D = 0, the apparent horizon and the spacetime singularity 
are formed on the same hypersurface x + Xq — —1, and now the singularity is marginally 
naked. In the last case, it can be shown that 

Xg^Xo^O, (£> = 0), (4.42) 

and the corresponding solution is given by 
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S = x^/'e-^/\ y = {D = 0). (4.43) 

3x 

Similar to the case po = and a = 1 discussed above in this section, these solution may 
also represent critical collapse. To have a definite answer to this problem, we need to study 
perturbations of the "critical" solution, which is out of the scope of this paper, and we hope 
to return to this problem in another occasion. 

On the other hand, the mass function defined by Eg. ( p. 91) now takes the form, 

m(t,r) = ^^^^e". (4.44) 
18 

Thus, on the apparent horizon x = Xg, we have 

MBH{t) = ^—-^ TAHit), (4.45) 

which shows that, as rAH{t) +oo, the total mass of black hole becomes infinitely large. 
Similar to the cases considered above, to have a black hole with finite mass, we can make a 
"surgery" to the spacetime. By this way, we can see that the resultant solution will represent 
a collapsing ball with a finite radius tq, and its contribution to the total mass of black hole 
is given by 

Mbh = ^-—^ ro, {D > 0), (4.46) 

which is finite and non-zero for any given non-zero Vq. 

C) Case 1 < a < 3: In this case the spacetime singularities and apparent horizon 
are still given by Eqs.( [4.3^ and ([4.35|) , respectively. It can be shown that the curve that 



represent the apparent horizon in the (t, x)-plane now is similar to that given in Fig. 5, that 
is, in the present case the singularity at x + xq = —2/(3 — a) is initially naked. As the fiuid 
is collapsing until the moment tc given by Eg. ( [4. 361) , the apparent horizon starts to form. 
The corresponding Penrose diagram is given by Fig. 6. 
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FIG. 12. The curves of the functions Yi(x) and Y2{x), defined by Eq.( [4.39| ) versus x. The 
equation Yi{x) = Y2{x) in general has three real roots, xg, xiq and xu, with xn > xiq > xg. 

D) Case a = 3: In this subcase, the corresponding physical quantities are given by 

3(x + Xq) — 2 



P 



9rft2(x + Xq) 



-9(-3tmx + xo)2/3-^' ^^-^^^ 
which show that the spacetime is singular at 

a)t = 0, b) x + xo = 0. (4.48) 

The location of the apparent horizon is given by 

[Ro.R^g''^ = O) , (4.49) 



' + ^4/3 



[x + Xq) 



(x + xo)^/3 

from which we can see that in the (t, x) -plane it is given by the curved given in Fig. 7. 
Consequently, the solutions in this case represent gravitational collapse that always forms 
black holes. 

E) Case a > 3: In this case, following the same routine given above, it is not difficult 
to show that the spacetime singularity at x + xq = 0, formed due to the collapse of the 
perfect fluid, is also covered by an apparent horizon. It can be shown that in the last two 
subcases the mass of such formed black holes is always finite and non-zero. 
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D. Case Po > ^ 



In this case, the corresponding solutions are given by 

1 



y[x) 



-- i tan 

3 r 



-{X + Xo) 



+ a 



S'(x) = e-°^/='cos2/3 



P 



{X + Xq) 



(4.50) 



where /3 is given as that in Eq.( [4.24| ). Then, the energy density of the fluid is given by 



/3(2a-3){(3po-l 


^tan 


f (a; + Xq) 


} 


3a2rft2 |tan 


^{x + xo) 





X < tan 



'/3 



X + Xq) 



B 



(4.51) 



where A and B are given by Eg. ( |4.26| ). From the above expression we can see that p is 
non-negative only in certain regions and the spacetime is singular on various hypersurfaces. 
The solutions in this case cannot be interpreted as representing gravitational collapse of 
perfect fluid. 



V. CONCLUDING REMARKS 

In Sec. Ill, we have studied the self-similar solutions of the zeroth kind, and found that 
some represent cosmological models and some represent gravitational collapse, while the 
others have no physical meanings. The ones that represent gravitational collapse are given 
by p = 0, i.e., dust fluid. These dust fluid solutions always collapse to form black holes with 
finite and non-zero mass. 

In Sec. IV, self-similar solutions of both the first and the second kinds have been studied. 
In particular, it has been found that the self-similar solutions of the first kind (a = 1) with 
Po = or = 1/3 may represent critical collapse but in the sense that now the "critical" 
solution separates solutions that form black holes to the solutions that form naked singulari- 
ties. In this case the formation of black holes also starts with a mass gap. To show explicitly 
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that these solutions indeed represent critical collapse, the analysis of spectrum of pertur- 
bations of these "critical" solutions is needed, which are currently under our investigation. 
The solutions with po = and a > 3/2, the ones with < po < 1/3 and ai < a < a2, and 
the ones with po — 1/3 and a > 3 also represent gravitational collapse, and the collapse 
always forms black holes with finite and non-zero mass. The solutions with pq — and 
< a < 1 and the ones with po = 1/3 and < a < 1 represent the formations of spacetime 
singularities that are covered by apparent horizons at the beginning of the collapse and late 
become naked, while the ones with po = and 1 < a < 3/2 and the ones with po — 1/3 
and 1 < a < 3 represent the formations of spacetime singularities that are naked at the 
beginning of the collapse and late become covered by apparent horizons. All the rest of 
the solutions can be either considered as representing cosmological models with an initial 
spacetime singularity or have no physical meanings. 

In review of all the above, one can see that the BC solutions with self-similarity of the 
zeroth and second kinds seem irrelevant to critical phenomena in gravitational collapse, and 
the only possible candidates for critical collapse are those solutions with self-similarity of 
the first (homothetic) kind with po = or po = 1/3, given in Sec.IV. 
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APPENDIX 



The metric for spacetimes with spherical symmetry can be cost in the general form, 

ds^ = rl {e2*(*'^)t/t2 - e2*(*'^)dr2 - r^S{t, rfdn''] , (A.l) 

where dQ^ = dO"^ + sin(6')^(i(y9^, and ri is a constant and has dimension of length, /. Then, 
it is easy to show that the coordinates {x'^} = {t, r, 9, ip}, the Christoffel symbols, V^^, 
the Riemann tensor, RJ^^xy ^^e Ricci tensor, i?^,^, and the Einstein tensor, G^i,, are all 
dimensionless, while the Ricci scalar, R, has the dimension of /~^, and the Kretschmann 
scalar, / = R^'^'^^R^^ux, has the dimension of 



For the metric ( [A.ID , we find that the non- vanishing Christoffel symbols are given by 



^ 00 




pO _ ^ pO _ 

01 ~ ^,r^ ^ 11 - 


= e2(*- 




pO 

22 


= r'^Se 


~2$r. pO _ ^2 


5sin2^ 




^ 00 


= e2(*- 


tL = ^ 

A'' 01 1 


J- 11 




22 


= -r5e-2* (rS^r + S) , 


33 ~ 


—rS sin^ 6e~^^ 


p2 

02 


S,t 
S ' 


■p2 '"'S'.r + S 
12 ~ c ' 

rb 


p2 _ 

33 ~ 


- — sin 6* cos 0, 


03 


S,t 


■p3 ''^^,r + 


p3 _ 
23 - 


cos 6 


S ' 


rS ' 


sin 6' ' 



[rS,r + S) 



while the non-vanishing components of the Einstein tensor are given by 



,-2'!' 



Gtt 
Gtr 

G 

Goo 



{e^* [2r^SS^rr + rS^r + QS) - 2rS {rS^,. + S) + - e^* 



^2^2 

-4 l''S,tr - irS,r + S) - S,t (r*,, - 1)] 

To 

— - {e^* 2r5 (r5 , + 5) $ , + r5 , (r5 , + 25) + 5^ - e 
r^D^ L L 

-rV^ [255ti + 5i(5i-25<l>,i)]}, 
r5e-2(*+*) {e^* [r (5$,,, + S,rr) + r5$,, - ^,,) 
+ (r5, + 5) - ^,,) + 25,] 
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2'!' 



(A.2) 



(A.3) 
(A.4) 

(A.5) 



(A.6) 



where ( = d{ )/dx°', ect. 

A. Solutions with Self-Similarity of the zeroth Kind 

To study solutions with kinematic self-similarity of the zeroth kind, let us introduce two 
new dimensionless variables, x and r, via the relations 



X = In(^) = — t + ln{r), t = t, 



(A.7) 



or inversely 



t = r, r = 



(A.8) 



Then, for any given function f{t, r) we find that 



1 



f.t f,T f ,xi f,r f,xi 

r 

1 1 

f.tr if,xx f,Tx) ; f,rr ~;T {f,xx /,a 

f,tt = f,TT — '^f,TX + / ,xx- 



(A.9) 



Gtr 



Substituting these expressions into Eqs.( [A.3| )-( [A.6| ), we find that 

-r^e^-" iS,r - S,^) [2S i^,^ - ^,x) + (5, - 5,)]} , 

^ [S,xx - S,rx + (5. + S) i^^r - + - 5,.) - l)] , 

rb 



{' 



,2<l> 



(A.IO) 
(All) 



25$,^ {S^^ + S) + {S^^ + 2S) + S^- e^* 
r e [2S {S^-TT ~ 2S + S^^x) 
+ (S, - ^,,) (S, - 25$,, - 5, + 25$,,)]} , (A.12) 
Gee = Se-2(*+'^) {e^* [^$,,, + + - m (5$,, + + S) + S,, - S^J 

- {S^>,r + - 5*,. - ($ ,r " " + • (A.13) 
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For the solutions with self-similarity of the zeroth kind, the metric coefficients $, and S 
are functions of x only, 



$(r,a;) = ^(r, x) = <l>(a;), S{r,x) = S{x). 



(A. 14) 



Then, all the derivatives of these functions with respect to r are zero, and Eqs.( |A.l"0[ ) 
( |A.13| ) reduce to. 



Gtt 



{e^* [2y,, - 2(1 + y)^,, + + 4y + 1 - S-^e^^ 
-rV*y(2v^,. + l/)}, 



Gt. = - [y,cc - (1 + y) -y)- y'^,x] 



(A.15) 
(A.16) 



-r^e^* (2y,, - 21/$,, + Sy^)} , (A. 17) 

Gee = S2e-^(*+'^) {e^* + y,, + - v];^, + ^) _ (i + ^) (v]/^, _ y)] 

^,xx + y,x + {^,x - ^,x) + y) + y^]}, (A.18) 



2 2* 

— r e 



where 



For a perfect fluid, the energy-momentum tensor (EMT) takes the form. 



(A.19) 



(A.20) 



where, when the fluid is co-moving with the frame of the coordinates, its four- velocity, u^, 
is given by 



(A.21) 



Then, from the 01-component of the Einstein field equations G^i, = T^jy, and Eq.( |A.l6| ) we 
find that 



y,, - (1 + y) (^,. -y)- y<l>,x = 0, (Goi = T, 



01 y 



(A.22) 
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while the other components yield 

-L^ {e-^^ (2y$,. + (1 + yf - S-'e 
V2*(2v|/,, + y)y}, (Goo = Too) 



[rir) 



P 



[Tir 



[rir) 

„2, 



1 + y) (2$,. + y + l)-S-V* 
[2(1 + y)^,, + y{y - 2)]} , (Gn = Tn) , 

e-'* [^,.. + (^,. - + 2t/ + 1) - t/]} , (G22 = T, 



(A.23) 



(A.24) 



(A.25) 



Note that in writing the above equations, Eg. ( [A. 221) was used. To have Eqs.( [A.2"^ and 
( |A.25| ) be consistent, we must have 



+ (^,x - - 1) - yiy - 1) = 0, 

+ $,x. ($,x. - - 2) - (1 + yy + S-^e^-" = 0, (a = 0). 



(A.26) 
(A.27) 



B. Solutions with Self-Similarity of the Second Kind 

To study solutions with kinematic self-similarity of the second kind, let us introduce 
other two new dimensionless variables, x and r, via the relations 



2; = In 



l/a 



In(-t) 



(A.28) 



or inversely 



t = -e-\ r = e^-"/", 



(A.29) 



where a is a dimensionless constant. Then, for any given function f{t,r) we find that 



f,t = 7 (a/,r + f,x) 

at 



f,r f,xi 



f,tr 



atr 

f,tt = (aV.rr + 2a/,^^ + /,^.^. + aV,r + a/,. 



if,XX ~\~ Oif,Tx) ; f,7 



if, XX f,x) 



(A.30) 
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Substituting these expressions into Eqs. (|A.3|) - (|A.6|) , we find that 



G 



tr 



Gr 



-2'f 



{a'e'" + S,, {S,, + AS) - 25'^,, (5,,. + S) + - e^* 



arts 



+a [5.,. - (5, + 5) - 5, ($,, - 1)]} 



-2$ 
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{a^e^* [25$,, (5 , + 5) + 5, (5, + 25) + 52 - e 



2* 



-^e^* [255,, + 5, (5, - 25$,, + 2a5)] 



■e^* [25 {aS^rr + 25,,) + 5, (5 , - 25$,,) 
+5, (a5,, - 2a5$,, + 5, - 25$,, + 2a5)]} 



t2 



Gee = 4 {"'e-'"" [5$,„ + 5,,, + ($,, - ^,,) (5$,, + 5, + 5) - 5$,, + 5, 
2 

_^e-2* [5^,„ + 5,„ - (5^,, + 5,,) ($,, - - a)] 



[5 (a^,,, + 2^,„) + a5,, 



25. 



- (5*,, + 5,,) (a$,, - a^,, + - - a) 
-($,,-^,,)(5^,, + 5,,)]}. 



(A.31) 



(A.32) 



(A.33) 



(A.34) 



For the solutions with self-similarity of the second kind, the metric coefficients are also 
functions of x only, but now with x being given by Eq. ( [A.28| ) . Then, setting all the derivatives 
with respect to r zero, Eqs. dOll ) - (|A:3^) reduce to 

1 



Gtt 



,2{*-*) 



2„2<I' 



22/,, + y{?>y + 4) + 1 - 2(1 + 1/)^,, - S~'e 



(2^,, + y)2/, 



Gtr = — [y,x + {l + y){y- ^,,) - 2/$,,] , 



G 



1 r 



rr ^2 



2„2<I' 



2(1 + 2/)$,, + (1 + 2/)^ -5-^6 
^e2(*-*) [2y,, + y{3y- 2$,, + 2a)] 



(A.35) 
(A.36) 

(A.37) 
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Gee = S'e-''' + y,^ + - + y) + {1 + y) {y - ^,,.)] 



-2$ 



2+2 



'^,xx + y,x + y - + y) - - «) 



(A.38) 



For a perfect fluid, the EMT is given by Eqs. ([A.20| ) and (|A.21|) . Similarly, from the 
01-component of the Einstein field equations we find that 



y,,-{l + y) (^,, -y)- y^,, = 0, (Gqi = T, 



01 > 



(A.39) 



— 



-2'!' 



Krf Lr"' 



2+2 



while the other components yield 

2^$,,. + (1 + yf - S-'e'"] } , (Goo = Too) 
2(l + y)$,. + (l + l/)2-S-V* 
2(1 + y)^,,. + y^ + 2{a - l)y] } , (Gn = Tn) , 
P = {;^e-^* + ($,. - VP,,. + 2y)] 

[^,,, + (^,, _ + 2y + a + 1) + (a - , (G22 = 



aH 



(A.40) 



(A.41) 



(A.42) 



where in writing the above equations, Eq. ( |A.31^ ) was used. 

When a 7^ 1, in the expressions of p the term that is proportional to r"^ has different 
power-dependence on r from the term that is proportional to when these expressions 
are written in terms of r and x, since 



(A.43) 



Then, the two expressions of Eqs. (|A.41|) and ( |A.42| ) are equal only when 



"if, XX + (^,x - ^,x) + («-!) (^,x -y)-y^ = 0, 
^,xx + ($,x - ^,x - 2) - (1 + yf + S-^e^-" = 0. 



(A.44) 
(A.45) 



When a = 1, all the terms in the expressions of p have the same power-dependence on 
r, and the two expressions of Eqs.( |A.4"l| ) and (|A.42|) are equal, provided that 
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+ - - 2) - (1 + vY + ^-'e 



,2(x+*-$) 



0, (a = l). 



(A.46) 



From the above equations we can see that a solution that satisfies Eqs.( |A.4"4D and (|A.45|) 
with a = 1 is also a solution of Eq. ( |A.4(^ ), but not the other way around, that is, a solution 
of Eq. (|A.46|) doesn't necessarily satisfy Eqs.( |A.4^ ) and (|A.45|) . 
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